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Abstract — Consensus  algorithms  have  generated  a  lot  of  interest 
due  to  their  ability  to  compute  globally  relevant  statistics  by  only 
exploiting  local  communications  among  sensors.  However,  when 
implemented  over  wireless  sensor  networks,  the  inherent  iterative 
nature  of  consensus  algorithms  may  cause  a  large  energy  con¬ 
sumption.  Hence,  to  make  consensus  algorithms  really  appealing 
in  sensor  networks,  it  is  necessary  to  minimize  the  energy  neces¬ 
sary  to  reach  a  consensus,  within  a  given  accuracy.  We  propose  a 
method  to  optimize  the  network  topology  and  the  power  allocation 
over  each  active  link  in  order  to  minimize  the  energy  consumption. 
We  consider  two  network  models:  a  deterministic  model,  where 
the  nodes  are  located  arbitrarily  but  their  positions  are  known, 
and  a  random  model,  where  the  network  topology  is  modeled  as 
a  random  geometric  graph  (RGG).  In  the  first  case,  we  show  how 
to  convert  the  topology  optimization  problem,  which  is  inherently 
combinatorial,  into  a  parametric  convex  problem,  solvable  with 
efficient  algorithms.  In  the  second  case,  we  optimize  the  power 
transmitted  by  each  node,  exploiting  the  asymptotic  distributions 
of  the  eigenvalues  of  the  adjacency  matrix  of  an  RGG.  We  further 
show  that  the  optimal  power  can  be  found  as  the  solution  of  a 
convex  problem.  The  theoretical  findings  are  corroborated  with 
extensive  simulation  results. 

Index  Terms — Consensus  algorithms,  minimum  energy  con¬ 
sumption,  random  graphs,  sensor  networks,  topology  control. 


I.  Introduction 

AVERAGE  consensus  algorithms  have  received  consider¬ 
able  attention  in  recent  years  because  of  their  ability  to 
enable  globally  optimal  decisions  using  only  local  exchange  of 
information  among  nearby  nodes  [1]— [3] .  The  price  paid  for 
this  simplicity  and  the  underlying  decentralized  philosophy  is 
that  consensus  algorithms  are  inherently  iterative.  As  a  con¬ 
sequence,  the  implementation  of  consensus  algorithms  over  a 
wireless  sensor  network  (WSN)  requires  an  iterated  exchange 
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of  data  among  the  nodes,  which  might  cause  an  excessive  en¬ 
ergy  consumption.  This  must  be  contrasted  with  a  centralized 
strategy  where  there  is  a  sink  node  that,  after  collecting  all  the 
observations  from  the  sensors  (perhaps  over  multiple  hops),  is 
virtually  able  to  compute  the  desired  statistic  in  a  single  shot. 
Hence,  to  make  consensus  algorithms  practically  appealing  in 
a  sensor  network  context,  it  is  necessary  to  minimize  the  en¬ 
ergy  consumption  necessary  to  reach  consensus.  Clearly,  the 
network  topology  plays  a  fundamental  role  in  determining  the 
convergence  rate  [4].  It  is  well  known  that,  as  the  network  con¬ 
nectivity  increases,  so  does  the  rate  of  convergence.  However, 
having  a  densely  connected  network  requires  a  high  power  con¬ 
sumption  to  guarantee  reliable  direct  links  between  many  nodes. 
In  principle,  having  a  fully  connected  network  is  equivalent  to 
having  as  many  sink  nodes  as  sensors,  so  that  the  convergence 
time  of  fully  connected  networks  is  minimum.  However,  the 
power  consumption  necessary  to  maintain  a  fully  connected  net¬ 
work  is  also  maximum.  On  the  other  hand,  a  minimally  con¬ 
nected  network  entails  low  power  consumption  to  maintain  a 
few  links,  but,  at  the  same  time,  it  requires  a  large  convergence 
time.  Since  what  really  matters  in  a  WSN  is  the  overall  power 
spent  to  achieve  consensus,  this  paper  addresses  the  problem 
of  finding  the  optimal  network  topology  that  minimizes  this 
overall  power  consumption,  taking  into  account  convergence 
time  and  transmit  power  necessary  to  establish  reliable  links 
jointly.  The  search  for  the  optimal  topology  is,  per  se ,  a  combi¬ 
natorial  problem  whose  solution  becomes  prohibitive  even  for 
small  scale  networks.  Nevertheless,  we  will  present  a  relaxation 
technique  that  converts  the  problem  into  a  convex  problem,  with 
minimum  performance  loss. 

If  the  links  among  the  nodes  are  symmetric,  or,  equivalently, 
if  the  graph  describing  the  network  topology  is  undirected,  the 
convergence  rate  can  be  lower  bounded  through  the  so  called 
algebraic  connectivity ,  defined  as  the  second  smallest  eigen¬ 
value  of  the  graph  Laplacian  [5].  For  this  reason,  there  have 
been  works  aimed  at  maximizing  the  algebraic  connectivity  of 
a  given  undirected  graph  by  a  suitable  choice  of  the  weights  as¬ 
sociated  to  each  edge  [6],  [7].  Alternatively,  in  [8]  it  was  shown 
how  some  network  topologies,  such  as  small  world  graphs,  can 
greatly  increase  the  convergence  rate.  In  [9],  the  authors  show 
that  nonbipartite  Ramanujan  graphs  constitute  a  class  of  topolo¬ 
gies  maximizing  the  convergence  rate.  In  [10]  it  was  shown  how 
to  add  edges  from  a  given  set  to  a  graph  in  order  to  maximize  its 
algebraic  connectivity.  Other  works,  for  example  [11],  proposed 
strategies  to  improve  the  convergence  rate  of  gossip  algorithms 
through  geographic  routing. 
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However,  in  all  these  works,  the  focus  is  always  on  conver¬ 
gence  time  and  there  is  no  cost  associated  with  the  establish¬ 
ment  of  the  graph  topology.  Conversely,  in  our  work,  since  the 
graph  represents  a  real  network,  we  consider  as  the  optimiza¬ 
tion  metric  the  power  consumption  necessary  to  maintain  re¬ 
liable  communication  links  among  the  nodes,  taking  into  ac¬ 
count  the  radio  propagation  model,  multiplied  by  the  number  of 
iterations  necessary  to  achieve  consensus.  From  this  perspec¬ 
tive,  enforcing  a  small  world,  scale-free,  or  Ramanujan  graph 
topology  may  not  be  the  best  choice  for  a  wireless  network, 
whose  topology  should  depend  on  the  application  at  hand  as 
well  as  on  the  propagation  model. 

In  the  wireless  communication  network  context,  there  have 
been  works  on  optimizing  the  network  topology  in  order  to  min¬ 
imize  the  power  consumption  necessary  to  guarantee  connec¬ 
tivity,  e.g.,  [12],  [13].  These  works  concentrate  on  the  communi¬ 
cation  task  and  do  not  make  any  specific  reference  to  the  running 
application.  However,  it  is  now  well  established,  in  the  sensor 
network  context,  that,  whenever  possible,  an  efficient  wireless 
sensor  network  design  should  take  into  account  the  specific  goal 
of  the  sensor  network  [14].  For  this  reason,  we  focus  on  the 
achievement  of  consensus  in  a  WSN.  It  is  well  known  that  one 
of  the  most  crucial  parameters  in  WSNs  is  energy  consumption, 
because  in  many  contexts  it  is  hard  to  recharge  the  batteries  or 
scavenge  energy  from  the  environment.  For  this  reason,  in  most 
applications,  minimizing  energy  consumption  is  more  appro¬ 
priate  than  minimizing  convergence  time  (although,  there  are 
important  applications  where  the  latter  could  be  more  impor¬ 
tant).  In  [15],  the  minimum  energy  consumption  problem  was 
studied,  assuming  a  common  transmit  power.  As  shown  in  [15], 
there  typically  is  an  optimum  power  that  minimizes  the  energy 
necessary  to  achieve  consensus  within  a  prescribed  accuracy. 

In  this  work,  we  generalize  the  initial  idea  suggested  in 
[15]  and  we  propose  a  method  for  optimizing  the  network 
topology  and  the  power  allocation  across  every  link  in  order 
to  minimize  the  energy  necessary  to  achieve  consensus.  We 
consider  two  classes  of  networks:  a)  deterministic  topologies, 
with  arbitrary,  but  known,  node  locations,  and  b)  random 
geometries,  with  unknown  node  locations,  modeled  as  random 
variables.  In  the  deterministic  case,  we  optimize  both  topology 
and  power  allocation.  Differently  from  [6],  we  do  not  assume 
any  prior  topology,  as  the  topology  comes  out  as  a  result  of  the 
optimization.  Topology  optimization  is,  in  general,  a  combina¬ 
torial  problem  and  hence  an  NP-hard  problem  (recall  that  an 
undirected  graph  composed  of  n  nodes  may  assume  2n(n-1)/2 
topologies).  To  tackle  this  issue,  we  propose  a  relaxation  tech¬ 
nique  that  allows  us  to  formulate  topology  optimization  as  a 
convex  parametric  problem.  Then,  we  show  that  the  effect  of 
this  relaxation  on  the  performance  is  negligible. 

In  the  random  topology  case,  where  the  intemode  distances 
are  unknown,  we  show  how  to  optimize  the  single  (common) 
transmit  power,  modeling  the  network  topology  as  a  random 
geometric  graph,  a  model  suitable  for  wireless  networks.  We 
provide  both  theoretical  and  simulation  results,  exploiting  the 
theory  of  the  eigenvalues  of  random  geometric  graphs. 

The  paper  is  organized  as  follows.  In  Section  II,  we  briefly  re¬ 
view  the  consensus  algorithm.  In  Section  III,  we  introduce  our 
communication  model  and  formulate  the  optimization  problem. 


Section  IV  is  devoted  to  topology  optimization  for  arbitrary 
networks.  In  Section  V,  we  start  by  providing  a  closed  form 
expression,  albeit  approximate,  for  the  algebraic  connectivity 
of  a  random  geometric  graph.  Then  we  use  this  expression  to 
proceed  with  the  topology  optimization  for  random  geometric 
graphs.  The  analytical  findings  are  corroborated  with  extensive 
simulation  results. 


II.  Brief  Review  of  Consensus  Algorithms 


Let  us  consider  a  wireless  network  composed  of  n  sensors. 
The  network  topology  can  be  represented  as  an  undirected  graph 
G  —  {V,  E}  where  V  denotes  the  set  of  n  vertices  (nodes) 
Vi  and  E  C  V  x  V  is  the  set  of  bidirectional  edges  (links) 
dj  =  eji  connecting  vt  and  Vj.  Furthermore,  let  A  be  the 
n  x  n-dimensional  symmetric  adjacency  matrix  of  the  graph  G , 
with  elements  =  1  if  aj  E  E  and  atJ  =  0  otherwise.  Ac¬ 
cording  to  this  notation  and  assuming  no  self-loops,  i.e.,  an  =  0, 
Vz  =  1, . . . ,  77,  the  degree  of  node  Vi  is  defined  as  deg(^)  = 

n  n 

aij  —  Y  aji •  The  degree  matrix  D  is  defined  as  the  di- 
j= i  j-i 

agonal  matrix  whose  diagonal  entry  is  da  =  deg(v4*).  Let  Aft 
denote  the  set  of  neighbors  of  node  z,  so  that  |A*|  =  deg(z^)1. 
The  Laplacian  matrix  of  the  graph  G  is  the  n  x  n  symmetric 
matrix  L  :=  D  —  A,  whose  entries  are 


_t  _  J  deg(vi)  if  j  =  i 
'  -  jl  -  l  -aij  if 


(1) 


Given  a  set  of  measurements  ^(0),  collected  by  node  z,  for 
z  =  1 , . . . ,  77,  the  goal  of  the  consensus  algorithm  is  to  allow 
every  node  to  compute  a  globally  optimal  function  of  the  mea¬ 
surements,  say  ,  through  a  decentralized  mechanism  that  does 
not  require  the  presence  of  a  sink  node,  but  builds  on  the  interac¬ 
tion  among  nearby  nodes  only.  Let  us  consider,  for  simplicity, 
the  case  where  x*  is  the  average  of  the  measurements.  In  this 
case,  reaching  an  average  consensus  can  be  seen  as  the  mini¬ 
mization  of  the  disagreement  between  the  states  Xi  of  the  inter¬ 
acting  nodes.  One  of  the  nice  properties  of  the  Laplacian  is  that 
the  disagreement  can  be  expressed  as  a  quadratic  form  built  on 
the  Laplacian  [16] 

j(x)  ;=  \  Y^^xi~x^2 

i= i  jeAfi 

^  n  n  ^ 

=  I  E  J2  M®*  -  X3?  =  o  xTLX-  (2) 

i=  1  j= 1 


An  important  property  of  the  Laplacian  is  that  it  has,  by  con¬ 
struction,  a  zero  eigenvalue,  whose  multiplicity  is  equal  to  the 
number  of  connected  components  of  a  graph.  Hence,  a  graph  is 
connected  if  its  zero  eigenvalue  has  multiplicity  one  [5].  Fur¬ 
thermore,  if  the  graph  is  connected,  the  eigenvector  associated 
with  the  null  eigenvalue  is  the  vector  1,  composed  of  all  ones. 

The  minimization  of  the  quadratic  form  in  (2)  can  be  achieved 
using  a  simple  steepest  descent  technique.  In  continuous  time, 
the  minimum  of  (2)  can  be  reached  through  the  following  up¬ 
dating  rule  [2] : 

x(t)  =  —Lx(t)  (3) 

|  •  |  we  denote  the  cardinality  of  the  set. 
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initialized  with  ar(0)  :=  xq2.  By  construction,  the  eigenvalues 
of  L  are  all  nonnegative.  Hence,  the  convergence  of  (3)  is  guar¬ 
anteed.  In  particular,  if  the  graph  is  connected,  the  state  vector 
x(t)  in  (3)  converges  to  the  projection  of  the  initial  state  xq  onto 
the  nullspace  of  L,  spanned  by  the  vector  1,  i.e.,  [2] 

x(t)  — y  -  1  lTx0.  (4) 

t— ►  oo  n 


This  corresponds  to  having  all  the  nodes  converging 
to  a  consensus  on  the  average  of  the  initial  observa¬ 
tions  Moreover,  the  convergence 

rate  is  lower  bounded  by  the  slowest  decaying  mode, 
i.e.,  by  the  second  smallest  eigenvalue  of  L,  \2{L), 
also  known  as  the  algebraic  connectivity  of  the  graph. 
More  specifically,  if  the  graph  is  connected,  the  dynamic 
system  (3)  converges  to  consensus  exponentially  [2],  i.e., 
|| x(t)  —  x*l\\  <  ||^(0)  —  x*l\\0(e~rt),  with  r  =  X2 (L). 
Defining  the  disagreement  vector  as  A x(t)  =  x(t)  —  x*l9  we 
can  write  ||Ax(t)||  <  || Ax(0)\\O(e~rt).  As  a  consequence, 
the  convergence  time  Tc  can  be  defined  (see  also  [2])  as  the 
time  necessary  for  the  slowest  mode  of  the  dynamical  system 
(3)  to  be  reduced  by  a  factor  7  «  1,  i.e.,  the  time  for  which 
||Aa:(Tc)||  <  7||Aar(0)||.  Hence,  we  can  set 


log(7) 

A 2(L) ' 


(5) 


III.  Optimization  Criterion 


It  is  well  known  that  one  of  the  most  critical  issues  in  WSN  is 
energy  consumption.  Neglecting  for  simplicity  the  power  spent 
for  processing  operations  with  respect  to  the  power  to  be  used 
to  enable  wireless  communications,  the  overall  power  spent  to 
reach  consensus  is  the  product  between  the  sum  of  the  power 
Ptot  necessary  to  establish  the  communication  links  among  the 
nodes  and  the  number  of  iterations  2Vit  necessary  to  achieve 
consensus.  The  exchange  of  information  among  the  nodes  is 
supposed  to  take  place  in  the  presence  of  a  slotted  system,  with  a 
medium  access  control  (MAC)  mechanism  that  prevents  packet 
collisions.  Denoting  by  Ts  the  duration  of  a  time  slot  unit,  the 
number  of  iterations  is  then  approximately3  N[t  =  Tc/Ts. 

Introducing  the  power  coefficients  ptj,  i  /  j,  denoting  the 
power  used  by  node  i  to  transmit  to  node  j  with  ptj  =  pji ,  and 
the  binary  coefficients  atj  assessing  the  presence  (atj  =  1)  of  a 
link  between  nodes  i  and  j  or  not  (a^  =  0),  the  power  spent  by 
the  whole  network  in  each  iteration  is  then  Ptot  —  Yi,  j  &ijPij  • 
Using  (5),  our  goal  is  to  minimize  the  following  performance 
metric: 


s  =  PtotVt  =  K 


En  \—\n 

i= 1  Yj= 1  aijPij 

w 


(6) 


2The  discrete-time  counterpart  of  (3)  is  [2] 


x[k  +  1]  =  x[k]  —  eLx[k]  Wx[k] 


where  e  is  a  parameter  chosen  so  as  to  ensure  that  no  eigenvalue  of  W  is  greater 
than  one  in  modulus. 

3 We  neglect  rounding  errors,  which  tend  to  vanish  if  the  duration  of  the  up¬ 
dating  time  slot  is  small  with  respect  to  the  convergence  time. 


where  K  incorporates  all  irrelevant  constants.  This  metric  is 
proportional  to  the  integral  of  power  consumption  with  respect 
to  time  and  then  it  represents  an  energy  consumption.  In  (6), 
we  have  made  explicit  the  dependence  of  the  Laplacian  L  on 
the  vector  a  =  A(:)  containing  all  binary  coefficients  atj,  since 
finding  these  coefficients  is  the  goal  of  our  optimization.  More 
specifically,  our  goal  is  to  find  the  set  of  active  links,  i.e.,  the 
nonzero  coefficients  a^,  and  the  powers  pLJ  that  minimize  (6), 
under  the  constraint  of  maintaining  network  connectivity,  which 
entails  A2 (L(a))  >  0.  The  problem  can  then  be  formulated  as 
follows: 


En  \—\n 

i=  1  Yj= 1  aijPij 

mm  - \  (T(  \\ - 

a,p  A  2  (//(a)) 

s.t.  6  <  A 2(L(a))  [P.0] 

O'ij  £  {0?  1} 

Pij  >  0  =  (7) 


where  e  is  an  arbitrarily  small  positive  constant  used  to  prevent 
the  algebraic  connectivity  from  going  to  zero,  which  would  cor¬ 
respond  to  a  disconnected  network  and  p  is  the  vector  with  en¬ 
tries  p^ . 

Since  the  topology  coefficients  are  binary  variables,  [P.0]  is 
a  combinatorial  problem,  with  complexity  increasing  with  the 
size  n  of  the  network  as  2n(n_1)/2 .  Hence,  its  solution,  for 
medium/large  scale  networks  is  prohibitive.  Our  objective  is  to 
modify  [P.0]  in  order  to  turn  it  into  a  convex  problem,  with  neg¬ 
ligible  performance  losses. 

A  first  simplification  comes  from  observing  that  the  coef¬ 
ficients  dij  and  p^  are  not  independent  of  each  other.  Their 
dependence  is  indeed  a  consequence  of  the  radio  propagation 
model.  In  this  work,  given  the  complexity  of  the  topology  op¬ 
timization  problem,  we  assume  a  fairly  simple  communication 
model.  We  state  that  there  is  a  link  between  nodes  i  and  j,  and 
then  atj  =  1,  if  the  signal-to-noise  ratio  SNRy  at  the  receiver 
node  j,  when  node  i  transmits,  exceeds  a  minimum  value  sm -in, 
i.e.,  SNRy  >  sm[n.  If  we  denote  by  pRj  the  power  received  by 
node  j  when  node  i  transmits,  and  by  <r2  the  noise  power,  as¬ 
sumed  for  simplicity  to  be  the  same  at  each  receiving  node,  we 
have  dij  =  1,  if  PRj  >  smin<Tn  Pmin,  or  otherwise  dtJ  =  0. 
Assuming  flat  fading  channel  modeling,  we  use  the  following 
propagation  model: 


PRj  = 


Pij 

1  +  (rij/ro)* 


(8) 


where  rzj  is  the  distance  between  nodes  i  and  j,  and  77  is  the 
path  loss  exponent.  The  parameter  77  plays  the  role  of  a  scaling 
factor  or  reference  distance,  and  typically  corresponds  to  the  so 
called  Fraunhofer  distance,  such  that,  if  rjj  ro,  the  receiver 
is  in  the  transmit  antenna  far-field,  where  the  received  power 
is  inversely  proportional  to  7*7.;  conversely,  if  r\j  <C  ro,  the 
receiver  is  in  the  transmit  antenna  near-field,  where  the  received 
power  is  approximately  equal  to  the  transmitted  one.  The  unity 
term  in  the  denominator  of  (8)  is  used  to  avoid  the  unrealistic 
situation  in  which  the  received  power  could  be  greater  than  the 
transmitted  one.  Given  the  propagation  model  (8),  the  relation 
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between  the  power  coefficients  p^  and  the  topology  coefficients 
a,ij  is  then 


1  if  Pij  ^  ^min 
0  otherwise 

In  the  following  sections  we  will  show  how  to  relax  this  rela¬ 
tion  in  order  to  simplify  the  solution  of  the  optimal  topology 
control  problem.  We  will  consider  two  scenarios:  a)  a  deter¬ 
ministic  topology,  with  arbitrary  node  geometry,  where  the  dis¬ 
tances  between  the  nodes  are  known,  and  b)  a  random  topology, 
where  the  positions  of  the  nodes  are  unknown  and  modeled 
as  random  variables.  In  case  a),  since  the  distances  among  the 
nodes  are  known,  we  optimize  the  power  allocated  to  each  link. 
This  strategy  is  also  instrumental  in  determining  the  topology, 
which  is  equivalent  to  finding  the  n(n  —  l)/2  entries  av;?-  of 
the  adjacency  matrix.  Conversely,  in  the  random  topology  case, 
since  the  distances  are  not  known,  we  assume  broadcast  com¬ 
munications,  and  look  for  the  optimization  of  the  coverage  ra¬ 
dius  of  each  node.  The  deterministic  topology  case  will  be  the 
subject  of  Section  IV,  whereas  the  random  case  will  be  studied 
in  Section  V. 


IV.  Optimal  Topology  and  Power  Allocation  for 
Arbitrary  Networks 


In  the  case  where  the  distances  between  the  nodes  are  known, 
the  optimization  criterion  amounts  to  solving  problem  [P.0], 
which  involves  a  combinatorial  strategy  that  makes  the  problem 
numerically  very  hard  to  solve,  especially  for  medium/large 
scale  networks.  The  relation  (9)  reduces  the  set  of  unknowns 
to  the  set  ptj ,  but  the  problem  still  retains  most  of  its  difficul¬ 
ties.  To  make  problem  [P.0]  to  be  numerically  tractable,  we 
introduce  a  first  relaxation  so  that,  instead  of  requiring  av;?-  to 
be  binary,  we  assume  azj  to  be  a  real  variable  belonging  to 
the  interval  [0,1].  Under  this  assumption,  problem  [P.0]  is 
not  combinatorial  anymore,  but  it  is  still  a  nontrivial  nonlinear 
constrained  problem.  The  first  important  contribution  of  this 
paper  is  to  propose  a  relaxation  technique  that  transforms  the 
previous  problem  into  a  convex  problem  that  can  be  solved  with 
well  established  and  efficient  numerical  tools.  We  achieve  this 
goal  by  first  introducing  the  following  relationship  between  the 
coefficients  atj  and  the  distances  r\j : 


1  +  (' rij/rCij)c 


(10) 


where  a  is  a  positive  coefficient  and  rCij  is  the  coverage  radius, 
which  depends  on  the  transmit  power.  According  to  (10),  atj 
is  close  to  one  when  node  j  is  within  the  coverage  radius  of 
node  i,  i.e.,  rvj  <C  rCij ,  whereas  avy  is  close  to  zero,  when 
r ij  rCij .  The  switching  from  zero  to  one  can  be  made  steeper 
by  increasing  the  value  of  a. 

Given  the  propagation  model  (8),  the  coverage  radius  rCij  is 
related  to  the  power  ptj  and  the  minimum  power  required  for 
reliable  communication  pm[n ,  as  follows: 


r 


to 


Pij 


(11) 


Plugging  (1 1)  in  (10),  the  coefficients  atj  can  be  written  explic¬ 
itly  in  terms  of  the  power  coefficients  as  follows: 


aij  ~  aij(Pij)  ~ 


rfi  {Pij  -Pmin)01^ 


ro(Pi 


)a/t]  i 

!/  1  zi-Cmin 


(12) 


This  relation  can  be  also  inverted  to  find  the  coefficients  ptj  as 
a  function  of  a^9  as  follows: 


/  a..  W" 

Pij  —  Q.ip'ij)  —  Pmin  T-  k\  Q,  -  ■  J  (1*^) 


rV . 

with  ki  =  PminAf .  Expression  (12)  becomes  our  relaxed  ver- 

ro 

sion  of  (9)  and  it  allows  us  to  reduce  the  set  of  variables  to  the 
only  power  vector  p.  Consequently,  problem  [P.0]  can  be  re¬ 
laxed  into  the  following  problem: 


min 

p 

s.t. 


pT  1 

WW) 

e  <  A 2(i(p)) 
Pminl  <  P 


[P-1] 


(14) 


where,  thanks  to  (12),  the  Laplacian  is  now  written  explicitly 
in  terms  of  the  power  coefficients  ptj .  In  principle,  the  last  in¬ 
equality  in  (14)  makes  any  link  feasible.  But  this  does  not  imply 
that  the  final  network  will  be  fully  connected,  because  some  co¬ 
efficients  a,ij  might  turn  out  to  be  equal  zero,  implying  that  the 
link  between  node  i  and  j  is  not  active.  The  first  important  re¬ 
sult,  related  to  the  solution  of  (14),  is  the  following. 

Theorem  1:  Given  the  propagation  model  in  (8),  using  the 
relations  (12)  between  the  topology  coefficients  and  the 
power  terms  ptJ ,  problem  [P.l]  is  a  convex-concave  fractional 
problem  if  p  >  a. 

Proof:  Let  us  consider  the  objective  function  in  [P.l]. 
The  numerator  of  [P.l]  is  clearly  a  convex  function  of  p.  We 
only  have  to  prove  the  concavity  of  the  algebraic  connectivity 
A2 {L{p))  with  respect  to  the  transmit  powers.  As  a  first  step, 
we  prove  that  (Hj{ptj)  is  a  concave  function  of  ptJ.  Then,  we 
use  this  to  show  that  A2(  ■ )  is  a  concave  function  of  p. 

(0  cbij{pij)  is  a  concave  function  of  ptj :  We  compute  the 
second-order  derivative  of  the  function  in  (12)  with  respect  to 
Pij 

-  =  M(a  -  v)h  ~  +  v) (Pij  ~  Pmm)a/ri] 

(15) 


where  the  constants  k2,  ki  and  are  given  by 


k2  =  Tq  ,  k3  = 


fc4  =  rfjp 


a/r] 
min  * 


ahhiPii  -ffmin)a/>?  2 

[k4  +  k2 (pij  -  Pmin)a/v]3  v2  ’ 


Note  that  ^2^4  >  0;  since  ptj  >  pm[n  Vz,  j,  we  also  have 
ks  >  0.  From  (15),  we  see  that,  if  a  <  rj,  the  second-order 
derivative  is  always  nonpositive  and  then  atj  (pzj )  is  a  concave 
function  of  ptj . 

(H)  ^2(*)  is  concave  in  p\  Exploiting  the  properties  of  the 
Laplacian  [16],  we  can  write  the  quadratic  form  associated  to 
L  as  in  (2).  From  (12),  we  also  note  that  each  coefficient  atj 


P  min 
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depends  only  on  the  corresponding  link  power  ptj9  and  not  on 
the  other  link  powers.  Given  any  pair  of  power  vectors  p W  and 
p(2\  let  p  =  /3p^  +  (1  —  ft)p(2\  with  0  <  /3  <  1.  Hence,  we 
have 


^  n  n 

xrL(a(p))x=-  J2  J2  ~xif 

i=l  j=l 
n  n 

=  2  ^  (^l'}  +  (X  “  “  xj)2 


i=  1  7=1 

n  n 


>- 

“2 


+  (Pi?)]  (xi~xj)2 

i  =  lj=l 

=  /3xtL  (V1^  x  +  (1—  (3)xTL  (^a  (p^2^jj  x 

(16) 


where  the  inequality  follows  from  the  concavity  of  atJ  (ptj ) . 

We  further  recall  that  the  algebraic  connectivity  A2  (L)  is  the 
solution  of  the  following  positive  semidefinite  program  (SDP) 

[5]: 


X2(L)  =  min  xTLx.  (17) 

®_L1,||*||=1 

Hence,  using  (16),  we  obtain 

A2  (l  (/?p(i)  +  (i-/3>(2))) 

>  A2  (j3L  (p(1))  +  (1  -  f3)L  (p(2))) 

>  px2  (L  (p(1)))  +  (1  -  P)X2  (L  (p(2))) 

where  the  last  equality  follows  from  a  majorization  theorem  for 
eigenvalues4.  Thus,  the  algebraic  connectivity  is  a  concave  func¬ 
tion  of  the  vector  of  transmit  powers.  Hence,  (14)  is  a  convex- 
concave  fractional  problem  since  it  is  the  ratio  of  a  convex  and 
a  concave  function  of  p.  ■ 

Since  (14)  is  a  convex-concave  fractional  problem,  we  can  use 
one  of  the  methods  that  solve  quasi-convex  optimization  prob¬ 
lems,  see,  e.g.,  [18],  [19].  For  example,  we  can  use  the  nonlinear 
parametric  formulation  proposed  in  [19].  To  do  so,  we  introduce 
the  following  function: 


with  x *  E  ©,  if,  and  only  if 

h*  =  h(p*,x*)  =  min {f{x)  —  p* g{pc)  :  x  €  ©}  =  0 

where  h(p*9x*)  means  that  for  p  =  p*  the  minimum  of 
{fix)  —  p*g(x)  :  x  E  ©}  is  taken  on  at  x  =  x*. 

Before  applying  Theorem  2,  it  is  useful  to  further  convert  the 
convex-concave  optimization  problem  [P.l]  into  the  following 
parametric  problem: 

min  pTl  -  p\2(L(p)) 
s.t.  6  <  A 2(L(p))  [P.2] 

lPmin<p.  (18) 

By  Theorem  1,  A 2(L(p))  is  a  concave  function  of  p.  Hence,  the 
objective  function  in  [P.2],  as  sum  of  convex  functions,  is  a 
convex  function.  The  constraint  sets  are  convex.  Then,  problem 
[P.2]  is  a  convex  parametric  problem,  whose  solution  is  a  func¬ 
tion  of  the  parameter  p  that  controls  the  tradeoff  between  the 
global  transmit  power  and  the  convergence  time.  Later  on,  we 
will  show  how  to  find  the  optimal  p. 

Since  problem  [P.2]  is  convex,  it  can  be  solved  using  nu¬ 
merically  efficient  convex  programming  tools.  However,  before 
applying  any  convex  tool,  it  is  worth  noticing  that  the  feasible 
set  in  [P.2]  is  not  compact.  Hence,  even  if  a  solution  exists,  in 
principle,  it  could  be  unreachable  in  finite  time.  To  overcome 
this  potential  drawback,  we  propose  next  an  alternative  formu¬ 
lation  of  [P.2]:  Instead  of  looking  for  the  set  of  power  coeffi¬ 
cients  pij,  and  then  for  the  avy ,  using  (13),  we  can  reformulate 
[P.2]  so  as  to  look  directly  for  the  variables  azj .  Then  the  opti¬ 
mization  problem  in  (18)  can  be  rewritten  in  terms  of  vector  a, 
as  follows: 


min  <f>(a)  —  /xA2(L(a)) 

a 

s.t.  e  <  A2(L(a))  [P-3] 

0  <  a  <  1  (19) 

where  </>(a)  =  Y^i-i  Y^j= 1  Q(aij )•  We  verify  next  that  this 

problem  is  still  convex.  To  study  the  behavior  of  <j>(a)9  we  com¬ 
pute  the  second-order  derivative  of  q{aij)9  obtaining 


h(p)  =  min{prl  -  pX  2(L(p))  :peT} 


where  p  is  a  real  positive  parameter  and  T  =  {p  :  p  > 
IPmin,  A2(L(p))  >  e}.  In  order  to  find  the  solution  of  this 
problem  we  could  use  the  following  result,  proved  in  [19]  and 
[20]. 

Theorem  2:  Let  f(x)  and  g(x)  be  continuous  real-valued 
functions  Vx  E  ©,  where  ©  is  a  nonempty  compact  subset  of 
Un  and  g(x)  >  0  Vx  E  ©.  Then 


p  = 


M) 

g(x*) 


=  mm 


I  fix)  . 
\g(x)  ' 


X  G  0 


9} 


4Let  A,  B  be  Hermitian  matrices,  and  let  C  —  A  +  B.  Let  their  eigen¬ 
values  be  sorted  in  nondecreasing  order.  Then  the  vector  of  eigenvalues  A  (C) 
majorizes  the  vector  A(A)  +  A (B).  Since  XpL)  —  0  for  the  graph  Laplacian, 
it  follows  that  A 2(A  +  B)  >  A 2(A)  +  A 2(B).  See  [17,  Theorem  4.3.27]. 


Pqjaij)  _  k  n  f  (Hj_ 

dafj  1  a  V 1  —  ctij 


T)l  Oi  —  2  ^ 

(1  —  Clij)4 


ia~1+2aij)  ’ 


We  infer  that 


(ppXHj) 

dah 


>  0  7} 


a  2aa,ij  0. 


We  note  that,  if  rj  >  a,  then  d  >  0  for  0  <  atj  <  1, 

so  that  <f>(a)9  as  a  sum  of  convex  functions,  is  convex5.  Finally, 
the  algebraic  connectivity  A 2(L(a))  is  a  concave  function  of 
a ,  as  can  be  proved  following  the  same  steps  as  in  Theorem 

5Note  that  </>(a)  is  the  sum  of  functions  of  the  single  variables  ai3 .  Hence  the 
convexity  of  < f>{a j  can  be  studied  looking  at  the  convexity  of  the  single  functions 

q(aZJ). 
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1.  Then,  the  optimization  problem  [P.3]  is  also  a  convex  para¬ 
metric  problem,  perfectly  equivalent  to  the  original  problem  in 
(18),  since  the  change  of  variables  ptJ  =  q(aij)  in  (13)  ensures 
a  one-to-one  mapping  q  :  R  — ►  U,  for  0  <  atJ  <  1,  with 
image  covering  the  problem  domain  in  (18)  (see  [21,  p.  130]). 
Hence,  assuming  0  <  a  <  1  —  e  ,  with  e  an  infinitesimal  posi¬ 
tive  constant  so  that  the  feasible  set  of  (19)  T  =  {a  :  0  <  a  < 
1  — e  ,  A2 (L(a))  >  e}  is  a  compact  convex  set  in  [Rn(n-1)/2,  we 
ensure  that  the  set  of  minima  of  (19)  is  nonempty  and  by  the  con¬ 
vexity  of  the  problem  we  can  deduce  that  all  local  minima  are 
also  global.  Note  that  since  the  optimization  problem  in  (19)  is 
convex  an  optimal  solution  can  be  found  via  efficient  numerical 
tools.  Furthermore,  using  Dinkelbach’s  algorithm  [19],  based 
on  Theorem  2,  we  are  also  able  to  find  the  optimal  parameter  f  t 
in  [P.3].  More  specifically,  the  Dinkelbach’s  algorithm,  applied 
to  our  problem,  proceeds  through  the  following  steps: 

1.  Set  i  =  1  and  let  a,i  be  a  feasible  point  of  T,  with  ^  = 

Hai) 

X  2(Zr(&;))’ 

2.  Set  p  =  pi  and  find  a^i  €  T  that  solves  the  minimization 
problem  in  [P.3]; 

3.  If  |  h(^ai+ 1)  |  =  |  (j)(ai+1)-fj,X2(L(ai+1))  |<  e",  withe" 
an  arbitrarily  small  positive  constant,  stop  and  take  ai+ 1  as 
the  optimal  link  coefficient  vector;  otherwise,  set  i  =  i  + 1 , 

=  iMk and  g0 10  slep  2 

Since  the  topology  coefficients  atJ  obtained  in  this  way  are 
real  variables  belonging  to  the  interval  [0, 1],  to  obtain  the  net¬ 
work  topology,  it  is  necessary  to  quantize  them  to  convert  them 
into  binary  values,  1  or  0,  indicating  the  presence  or  absence 
of  a  link.  This  quantization  is  achieved  by  comparing  each  atj 
with  a  threshold  ath-  Of  course,  the  final  topology  will  depend 
on  the  threshold  value.  Moreover,  the  thresholding  operation 
will  also  affect  the  final  result  in  terms  of  convergence  time 
and  energy  consumption.  It  is  then  of  interest  to  check  how 
sensitive  the  final  topology,  as  well  as  convergence  time  and 
energy  consumption,  are  to  the  choice  of  ath-  In  the  ensuing 
section,  we  present  some  numerical  results  to  shed  light  on  the 
resulting  topologies  and  their  dependence  on  the  propagation 
model  parameters. 

A.  Numerical  Examples 

Since  our  optimization  procedure  is  based  on  a  relaxation 
technique,  the  first  important  step  is  to  evaluate  the  impact  of 
the  relaxation  on  the  final  topology  and  performance. 

Example  1:  Comparison  Between  Exhaustive  Search  and  Re¬ 
laxed  Technique:  We  compare  now  the  topology  obtained  as  a 
solution  of  the  relaxed  problem  [P.3]  with  the  optimal  graph 
obtained  by  solving  directly  problem  [P.0]  using  an  exhaustive 
search  over  all  possible  topologies.  For  complexity  reasons,  of 
course  we  can  only  perform  this  comparison  for  small  scale  net¬ 
works.  We  consider  networks  of  n  =  4  and  6  nodes. 

To  provide  results  not  conditioned  to  a  specific  geographic 
node  deployment,  we  averaged  the  results  over  100  statistically 
independent  realizations  of  the  nodes  locations.  In  each  iteration 
we  compute  the  minimum  energy  Eop  reached  using  the  optimal 
exhaustive  search  over  all  possible  topologies,  and  the  energy 
Er  corresponding  to  the  network  topology  whose  coefficients 
dij  are  obtained  by  solving  problem  [P.3]  and  thresholding  the 


Fig.  1.  Optimal  average  energy  obtained  with  exhaustive  search  and  average 
energy  obtained  by  solving  the  relaxed  problem  [P.3]  versus  77,  for  different  n 
values. 


... 


Fig.  2.  Ratio  between  the  optimal  energy  obtained  with  exhaustive  search  and 
the  energy  obtained  by  solving  the  proposed  relaxed  problem  [P .  3] ,  versus  at h , 
for  different  values  of  n . 


result  with  the  threshold  ath  =  2  x  10-3.  In  Fig.  1,  we  report  the 
average  energies  Eop  and  Er  versus  77.  We  can  note  from  Fig.  1 
that  the  loss  in  terms  of  optimal  energy  due  to  the  relaxation  of 
the  original  problem  is  negligible  (the  energy  loss  is  less  than 
1.7  dB). 

Example  2:  Impact  of  Thresholding  Operation  on  Final 
Topology:  Clearly,  the  selection  of  the  threshold  ath  plays  a 
role  in  the  identification  of  the  final  topology.  To  evaluate  the 
impact  of  ath  on  the  final  topology,  in  Fig.  2  we  report  the 
ratio  versus  at h  for  different  numbers  of  nodes  n  and  for 

£-op 

7}  =  5.  The  energies  are  averaged  over  100  independent  node 
realizations.  We  can  observe  from  Fig.  2  that  there  is  a  wide 
range  of  values  of  ath  such  that  the  energy  loss  is  practically 
independent  of  ath-  This  shows  that  our  proposed  procedure  is 
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rather  robust  with  respect  to  the  choice  of  at h-  It  is  now  inter¬ 
esting  to  check  the  effect  of  the  threshold  at h  on  the  topology, 
considering  larger  scale  networks  (in  this  case,  we  can  only  run 
our  algorithm,  as  the  exhaustive  search  is  not  feasible).  As  an 
example,  in  Fig.  3  we  show  the  topologies  obtained  by  solving 
problem  [P.3],  for  a  network  composed  of  n  =  30  nodes,  using 
different  values  of  ath,  for  rj  =  6.  Comparing  the  four  cases 
reported  in  Fig.  3,  we  notice  that,  only  for  very  low  values  of 
the  threshold  [i.e.,  case  (d)],  we  appreciate  a  sensitive  change 
of  topology,  whereas  for  a  large  range  of  values  of  ath,  the  final 
topology  is  practically  the  same.  This  means  that  the  proposed 
algorithm,  in  spite  of  the  relaxation  step  and  the  subsequent 
quantization,  yields  rather  stable  solutions. 

The  previous  results  pertain  to  a  specific  realization  of  the 
node  locations.  To  provide  results  of  more  general  validity,  in 
Fig.  4,  we  report  the  average  value  of:  a)  fraction  of  active  links 

n(n- 1)'1 1  ’  b)  ^2  (L) ;  and  c)  the  average  energy  £r ,  as  a  function 
of  tne  threshold  ath-  The  averages  are  carried  out  over  100  inde¬ 
pendent  realizations  of  the  nodes  location.  From  Fig.  4,  we  ob¬ 
serve  that  there  is  an  interval  of  values  of  at h  (roughly,  between 
10“ 5  to  10“ 2)  for  which  we  obtain  a  strong  reduction  in  the 
fraction  of  active  links,  with  respect  to  the  situation  where  there 
is  no  threshold,  still  achieving  nearly  the  same  performance,  in 
terms  of  algebraic  connectivity  and  energy  consumption.  This 
is  indeed  an  important  result,  as  it  shows  that  the  relaxed  algo¬ 
rithm  is  weakly  sensitive  to  the  choice  of  at h- 


Fig.  4.  Average  value  of  (a)  fraction  of  active  links,  (b)  A 3(L).  (c)  energy 
versus  threshold  value  for  rf  —  6 . 


Example  3:  Impact  of  Propagation  Parameters  on  Final 
Topology :  It  is  also  interesting  to  look  at  the  change  in 
topology  as  a  function  of  the  radio  communication  model.  To 
this  end,  in  Fig.  5,  we  plot  the  optimal  topologies  achieved  for 
the  same  node  locations  as  in  Fig.  3,  but  pertaining  to  different 
path  loss  exponents  77,  for  a  given  threshold.  Interestingly,  we 
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notice,  as  expected,  that,  as  the  propagation  loss  increases  (i.e., 
r )  increases),  the  network  tends  to  become  more  and  more 
sparse. 

The  results  shown  in  Fig.  5  refer  to  a  single  realization  of  the 
node  positions.  To  draw  conclusions  of  more  general  validity, 
we  averaged  over  100  statistically  independent  realizations  of 
the  node  locations,  for  networks  of  30  nodes.  In  Fig.  6,  we  report 
the  same  performance  metrics  as  in  Fig.  4,  but  now  as  a  function 
of  the  path  loss  exponent  77,  setting  ath  =  0.09.  From  Fig.  6, 
we  observe  that  when  the  attenuation  is  high  (i.e.,  77  is  large), 
reducing  the  number  of  links  (making  the  topology  sparser)  is 
more  important  than  reducing  convergence  time.  Conversely, 
when  the  attenuation  is  low  (i.e.,  77  is  small),  increasing  network 
connectivity  is  more  important  than  reducing  power  consump¬ 
tion.  This  behavior  sounds  reasonable  and  in  agreement  with 
intuition. 

V.  Optimal  Topology  for  Random  Geometric  Networks 

In  this  section,  we  remove  the  assumption  that  the  node  lo¬ 
cations  are  known  a  priori  and  model  the  network  as  a  random 
geometric  graph  (RGG).  In  such  a  case,  the  graph  connectivity 
properties  and  the  convergence  time  can  only  be  established  in 
a  probabilistic  sense,  asymptotically,  as  the  number  of  nodes 
tends  to  infinity.  We  refer  to  [22]  for  the  first  basic  result  about 
the  convergence  of  consensus  algorithms  over  random  graphs 
and  to  [23]  for  a  more  recent  generalization  of  the  convergence 
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Fig.  6.  Average  value  of  (a)  fraction  of  active  links,  (b)  A2(X).  (c)  energy 
versus  path  loss  77  for  ath  =  0.09. 

conditions  under  stochastic  disturbances.  As  shown  in  [22], 
the  rate  of  convergence  of  consensus  algorithms  in  the  random 
graph  case  is  dictated  by  the  expected  value  E[e~2TsX 2^]. 
In  the  following,  exploiting  the  concentration  properties  of  the 
eigenvalues  of  RGGs  [24],  [25],  we  will  show  how  to  relate 
the  convergence  time  to  the  expected  value  of  A2  (L).  This  link 
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will  be  fundamental  to  derive  the  optimal  coverage  radius,  and 
then  transmit  power,  that  minimize  the  energy  consumption 
necessary  to  achieve  consensus  over  RGGs.  To  achieve  this 
goal,  it  is  fundamental  to  recall  and  extend  some  results  about 
the  spectrum  of  RGGs. 


A.  RGGs 


A  random  graph  is  obtained  by  distributing  n  points  randomly 
over  the  d-dimensional  space  Rd  and  connecting  the  nodes  ac¬ 
cording  to  a  given  rule.  Let  Vn  =  {xi,x2, . . .  ,xn}  be  a  set 
of  d-dimensional  vectors  Xi ,  belonging  to  a  compact  set  ft  = 
[— 1/2, 1/2] d  of  Rd,  whose  entries  denote  the  coordinates  of  the 
nodes  in  Rd ,  and  let  ||  •  ||  be  the  Euclidean  norm  on  Rd.  The 
graph  topology  is  captured  by  the  adjacency  matrix  A  which, 
in  this  case,  is  a  random  matrix.  An  important  class  of  random 
matrices,  encompassing  the  adjacency  matrix  of  our  problem, 
is  the  so-called  Euclidean  Random  Matrix  (ERM)  class,  intro¬ 
duced  in  [26] .  Given  a  set  of  n  points  located  at  positions  x^  i  = 
l,...,n,annx?i  adjacency  matrix  A  is  an  ERM  if  its  generic 
(z,  j)  entry  depends  only  on  the  difference  Xi  —  Xj ,  i.e.,  = 

F(x{  —  Xj),  where  F  is  a  measurable  mapping  from  Rd  to  R. 
An  important  subclass  of  ERM  is  given  by  the  adjacency  ma¬ 
trices  of  the  so-called  RGGs.  In  such  a  case,  the  entries  atj  of 
the  adjacency  matrix  are  either  zero  or  one  depending  only  on 
the  distance  between  nodes  i  and  j,  i.e. 

(20 

1 0  otherwise 


where  r  is  the  coverage  radius.  This  is  a  particular  case  of  an 
ERM,  corresponding  to  having 


F(xi-Xj)  =  |  J 


if  || Xi  —  Xj ||  <  r 
otherwise 


(21) 


We  will  use  the  symbol  G(n,  r)  to  indicate  an  RGG  composed 
of  n  points,  with  coverage  radius  r. 

The  RGG  model  is  the  most  appropriate  to  capture  the 
topology  of  a  wireless  network,  as  it  basically  states  that  there 
is  a  link  between  two  nodes  only  if  they  are  within  the  coverage 
radius  of  each  other6. 

Next,  we  recall  some  of  the  most  important  properties  of 
RGG’s,  in  terms  of  connectivity  and  spectrum,  as  they  are  rele¬ 
vant  to  our  optimization  problem. 

1 )  Connectivity:  Some  interesting  results  on  the  asymptotic 
connectivity  of  random  geometric  graphs  have  been  derived  in  a 
seminal  work  by  Gupta  and  Kumar  [27]  who  proved  that,  given 
a  set  of  n  points  uniformly  distributed  within  a  unit  square  (i.e., 
d  =  2),  the  graph  is  connected  almost  surely  if  the  coverage 
radius  behaves  as 


r0(n) 


log  (n)  +  c„ 


7 m 


with  cn  — >  oo,  as  n  — >  oo.  Conversely,  if  cn  — ►  —  oo,  the  graph 
is  disconnected  almost  surely.  This  means  that  the  expression 

lo^n^  represents  a  threshold  distance.  In  the  following,  we 

6In  practice,  a  wireless  channel  is  also  affected  by  fading,  multipath,  and  shad¬ 
owing.  Hence,  the  presence  of  a  link  between  two  nodes  depends  on  two  sources 
of  randomness:  distance  and  fading.  In  this  work,  we  concentrate  on  the  single 
source  of  randomness,  given  by  the  distance  between  the  nodes,  but  considering 
both  sources  of  randomness  will  be  an  interesting  extension  of  this  paper. 


will  often  use  the  shorthand  notation  r+(n)  to  indicate  the  law 
^/lQg(n)+^  wjtj1  Cn  _ oo,  as  n  — ►  oo.  Hence,  a  coverage 
radius  ro(n)  behaving  as  r+(n)  represents  a  law  that  ensures 
connectivity  with  high  probability,  as  n  — ►  oo.  We  will  use  the 
notation  ro(n)  ~  r+(n)  to  indicate  such  a  behavior. 

In  [4],  [28],  it  has  been  shown  that  the  degree  of  a  RGG 
G(n,r)  of  points  uniformly  distributed  over  a  two-dimensional 
unit  torus7  is  equal  to 

d(n)  =  7r  r2n  (22) 


with  high  probability,  i.e.,  with  probability  1  —  1  /n2,  if  the  radius 
behaves  as  ro(n)  ~  r+(n).  This  implies  that  an  RGG  tends  to 
behave,  asymptotically,  as  a  regular  graph,  if  the  coverage  radius 
is  chosen  so  as  to  guarantee  connectivity  with  high  probability. 

We  are  primarily  interested  in  the  second  eigenvalue  of  the 
Laplacian,  L  =  D  —  A,  where  D  is  the  degree  matrix  and  A  is 
the  adjacency  matrix  [see  (1)].  From  (22),  D  =  7 rr2n/,  so  that 
we  only  need  to  investigate  the  second  largest  eigenvalue  of  A. 
Hence,  in  the  ensuing  section,  we  study  the  spectrum  of  A. 

2)  Spectrum  of  a  Random  Geometric  Graph:  In  [24],  [25], 
it  is  shown  that  the  eigenvalues  of  the  adjacency  matrix,  or  of 
the  transition  probability  matrix  8,  tend  to  be  concentrated,  as 
the  number  of  nodes  tend  to  infinity.  In  particular,  in  [24]  it  is 
shown  that  the  eigenvalues  of  the  normalized  adjacency  matrix 
An  =  A/n  of  an  RGG  G(n,  r),  composed  of  points  uniformly 
distributed  over  a  unit  bidimensional  torus,  tend  to  the  Fourier 
series  coefficients  of  the  function  F  defined  in  (21) 

F(z)  =  /  exp  (—2irjzTx)dx  (23) 

Jnr 


almost  surely,  for  all  z  =  [zi,z2\  E  Z2,  where  ftr  =  {x  = 
[xi,x2]t  E  R2  :  ||x||  <  r}.  Using  polar  coordinates,  i.e., 
x\  =  p  sin  0  and  x2  =  p  cos  0 ,  with  0  <  p  <  r  and  0  <  9  <  2ir, 
we  obtain 

f-r  t 

F(z)  =  /  /  exp  (— 2irjp(zi  sin  0  +  z2  cos  0))pdpd0. 

Jo  Jo 


This  integral  can  be  computed  in  closed  form.  Setting  zi  = 
A  sin  (j>  and  z2  =  A  cos  </>,  we  have 

F(A,(j>)=  /  /  exp  (— 27rjpAcos(t;))pdpdt; 

Jo  J-o 


with  £  =  0  —  cj).  Furthermore,  using  the  integral  expression 
for  the  Bessel  function  of  the  first  kind  of  order  k,  Jk(x)  = 

h  If  exP  0X  sin(0  -  we  §et 

pr 

F(A,(j>)  =  F(A)  —  2i r  Jo(27rpA)pdp. 

Jo 

Finally,  using  the  identity  JQn  v  Jo  (v)dv  =  uJ\  (u),  we  can  make 
explicit  the  dependence  of  F(A)  on  the  index  pair  [zi,z2] 

F(zi,z2)  =  2  J\  (2w rfffzA  ■  (24) 

V  Z1  +  z2  X  / 

7 A  torus  geometry  is  typically  used  to  get  rid  of  border  effects. 

8The  transition  probability  matrix  is  the  adjacency  matrix,  normalized  with 
respect  to  the  node  degree,  so  that  the  Ah  row  of  the  adjacency  matrix  is  divided 
by  the  degree  of  node  i. 
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Fig.  7.  Verification  of  the  inequalities  in  (26). 


This  formula  allows  us  to  rank  the  eigenvalues  of  An  =  A/ n. 
In  particular,  we  are  interested  in  the  second  largest  eigen¬ 
value  of  An.  Considering  that  the  minimum  coverage  radius 

ensuring  connectivity  behaves  as  r(n)  ~  i.e.,  it 

is  a  vanishing  function  of  n ,  we  can  use  the  Taylor  series 
expansion  of  F(zi1z2),  for  small  r.  Recalling  that,  for  small 
x,  J\(x)  =  x/2  —  x3/16  +  o{x 5),  we  can  approximate  the 
eigenvalues  as 


F(zi,z2)  =  7rr2 


ir^zl  +  zDr4 


+  o(r6).  (25) 


,\  =  Ittii 
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Fig.  8.  Asymptotic  analysis  and  simulation  results  for  some  eigenvalues  of  An 
and  for  X2(Ln)  versus  the  transmission  radius. 


shown  in  (22),  the  second  largest  eigenvalue  of  Ln  tends  asymp¬ 
totically  to 


M{Ln)  =  7T r2  -  C2(Ai).  (29) 

Thus,  the  algebraic  connectivity  of  the  graph  can  be  approxi¬ 
mated,  asymptotically,  as 


This  expansion  shows  that,  at  least  for  small  r,  the  largest  eigen¬ 
value  equals  i rr2  and  occurs  at  z±  =  Z2  =  0,  whereas  the  second 
largest  eigenvalue  corresponds  to  the  cases  {z\  —  1,^2  =  0) 
and  ( z\  —  0,  ^2  =  1).  More  generally,  we  can  check  numeri¬ 
cally  that,  for  r  <  1/2  and  A  >  1,  the  following  inequalities 
hold  true: 

7 rr2  >  rJi(27rr)  >  \  Ji(27rrA)|.  (26) 

A 

The  validity  of  these  inequalities  can  be  verified  from  Fig.  7, 
which  shows  the  three  terms  in  (26)  as  a  function  of  r  and  A. 

In  summary,  denoting  the  spectral  radius  of  An  as  Ci(^-n)  = 
maxi<^<n  ,  where  (A^n)}-^  is  the  set  of  eigenvalues 
of  A,  it  follows  that 

lim  Ci(^-n)  =  max  I  F(z)  1=  F( 0,0)  =  nr2  (27) 

n^-oo  z£  Z2 

while  the  second  largest  eigenvalue  of  An,  ( 2(An ),  converges 
to 

lim  C2(An)  =  F(1,0)  =  F(0,1)  =  rJi(27rr).  (28) 

n— »•  oo 

We  are  now  able  to  derive  the  asymptotic  expression  for  the 
second  largest  eigenvalue  of  the  normalized  Laplacian  Ln  = 
Dn  —  An ,  where  Dn  :=  D/n  is  the  normalized  degree  matrix. 
Because  of  the  asymptotic  property  of  the  degree  of  an  RGG, 


\2(L)  =  irnr2  —  nr  Ji(27rr).  (30) 

Since  the  previous  expressions  have  been  derived  in  the 
asymptotic  regime,  in  Fig.  8  we  compare  the  analytic  formulas 
for  the  first  and  the  second  largest  eigenvalues  of  An ,  as  given 
in  (27)  and  (28),  with  the  numerical  results  obtained  by  aver¬ 
aging  over  100  independent  realizations  of  RGG’s  composed 
of  n  =  1024  nodes.  We  can  notice  the  good  agreement  in 
Fig.  8  between  the  theoretical  expression  for  the  algebraic 
connectivity  A2 (Ln),  given  in  (29),  and  the  simulation  results. 

B.  Minimization  of  the  Energy  Consumption:  An  Analytic 
Approach 

We  can  now  exploit  the  previous  analytic  expressions  to  study 
the  energy  minimization  problem  for  RGG’s.  In  the  random 
topology  case,  since  the  distances  are  unknown,  we  cannot  op¬ 
timize  the  power  associated  with  each  link.  However,  we  can 
seek  the  common  transmit  power  that  minimizes  energy  con¬ 
sumption.  Thus,  in  the  random  setting  we  assume  a  broadcast 
communication  model,  where  each  node  broadcasts  the  value 
to  be  shared  with  its  neighbors.  In  the  lack  of  any  information 
about  distances  among  the  nodes,  we  assume  that  each  node  uses 
the  same  transmit  power.  In  this  case,  the  network  topology  can 
be  modeled  as  a  random  graph  model.  It  can  be  shown  [22], 
[23]  that  the  system  in  (3)  converges  to  consensus  almost  surely, 
i.e.,  Pr{  lim  x(t)  =  x*l}  =  1  assuming  that  each  node  has  a 

t — >oc 

coverage  radius  so  that  the  network  is  asymptotically  connected 


32 


with  probability  one  [14].  Then  the  rate  of  convergence  of  the 
dynamical  system  in  (3)  is  given  [22],  [23]  by  E[e~2TsX 2^)]. 
Note  that  defining  /( A2)  =  e~2TsX and  using  the  Taylor 
series  expansion  of  /( A2)9  at  the  point  A2  =  E[X 2]  =  m\2,  we 
can  write 


/(A2)  =  f(m\2)  +  /(1)(toa2)(A2  -  m\2) 
/(2)(m  A2)/\  _  n2 


+ 


-(A2  -  toa2) 


V2'  f(k\m a2)  fe 

+  2_^  fcj  (a2-toa2)  • 


k=3 


(31) 


so  that  g(r)  is  a  convex  increasing  positive  function  of  r.  Let  us 
now  study  the  behavior  of  A2(r).  Using 

—  =  rmJm-i(r)  for  m  =  0,1,2,... 

we  obtain 

A2(r)  =  r ^  —  27rrn[l  —  Jo(27rr)] 

with  A2(r)  >  0  for  ro(n)  <  r  <  0.5.  Moreover,  since 

dJm(r)  _m  _nio 

—  Jmy)  y  )  for  Tfl  —  0,1,2,... 

dr  r 


Consequently,  taking  the  expected  value,  we  get 

E[f(  A2)]  =  f(m  x2 )  +  - — y^~(t12  +  •  •  •  (32) 

denoting  with  cr\2  =  E[( X2  —  m\2)2]  the  variance  of  A2.  But 
since  all  central  moments  of  order  greater  than  one  of  the  eigen¬ 
values  tend  to  zero,  because  of  the  concentration  property,  we 
can  use  the  approximation 

E[ e~2T‘X2]  «  e~2T‘E[X2].  (33) 


As  a  consequence,  the  energy  spent  to  achieve  consensus  can 
now  be  approximated  as 


£  =  K _ — _ 

2E[X2(L(p))Y 


(34) 


This  is  the  performance  metric  we  wish  to  minimize  in  the 
random  scenario,  with  respect  to  the  single  unknown  p. 

In  particular,  using  the  asymptotic  expression  (30)  for  the  al¬ 
gebraic  connectivity,  we  can  introduce  the  following  metric 

F(r\  =  »Pmin[l  +  (r/r  0)^] 

[  l  nirr2  -  rnJ^wr)  '  V  ’ 


We  now  check,  numerically,  that  the  function  £(r)  given  in  (35) 
is  a  convex  function  of  r,  for  ro(n)  <  r  <  0.5,  where  ro(n)  ~ 
r+(n),  to  ensure  connectivity. 

Let  us  rewrite  (35)  as 


£{r) 


9_iT) 

A2(r) 


with  g(r)  =  npm[n[l  +  ( r/ro )T?]  and  A2(r)  as  in  (30).  The  first- 
and  second-order  derivatives  of  g(r)  are,  respectively 

/  dg(r )  r 17-1 

9  (r)  =  —7—  =  npminr) — —  >  0 
dr  r o1! 

d^n(v}  r7!  2 

9  (r)  =  ^  =  npminrj(V  -  1)— >  0 

9For  simplicity  we  drop  in  A2  the  dependence  on  L. 


we  have 

A 2(r)  =  ^  =  2tt^[1  —  Jo(27rr)]  +  47r2rnJi(27rr). 

(36) 


Observe  that  the  first  term  on  the  right-hand  side  (RHS)  of  (36) 
is  always  positive.  Furthermore,  if  r  <  1/2,  the  second  term  is 
also  positive,  since  Ji(27rr)  >  0  for  r  <  1/2.  Hence,  A2  (r)  >  0 
and  we  can  conclude  that  the  algebraic  connectivity  is  an  in¬ 
creasing  and  convex  function  of  r  for  ro(n)  <  r  <  0.5,  where 
r0(n)  ~  r+(n). 

We  can  now  compute  the  first-  and  second-order  derivatives 
of  the  energy  function  with  respect  to  r.  We  get 


e'(r)  = 


9  ir) 

A  2(0 


A2(r)g(r) 

A  Hr) 


and,  substituting  the  corresponding  expressions,  the  extremal 
points  can  be  obtained  by  solving  the  following  nonlinear 
equation: 


£  (r)  =  ai (r)  +  a2(r)J0(27rr)  -  a3(r)  Ji(27rr)  =  0 


with  ai(r)  =  7r[(r}  —  2)r^  —  2 rg],  a2(r)  =  27 r(rr?  +  ro*7)  and 
a3(r)  =  pr71-1 .  Furthermore,  the  second  derivative  of  £{r)  is 
given  by  the  equation  shown  at  the  bottom  of  the  page.  In  Fig.  9, 
we  report  £  (r)  as  a  function  of  r,  for  different  values  of  p. 
From  Fig.  9,  we  can  check  that  the  second  derivative  is  always 
positive  in  the  range  of  interest.  This  verifies  that  £{r)  is  in¬ 
deed  a  convex  function  of  r  defined  on  the  compact  convex  set 
r0(n)  <  r  <  0.5.  As  a  consequence,  we  can  state  that  there 
is  always  at  least  a  radius  r  that  globally  minimizes  the  energy 
consumption  in  a  RGG. 

Numerical  examples.  In  Fig.  10,  we  compare  the  value  of  £{r) 
obtained  by  our  theoretical  approach  and  by  simulation,  for  var¬ 
ious  values  of  the  path  loss  exponent,  2  <  p  <  6.  The  results 
are  averaged  over  100  independent  realizations  of  random  geo¬ 
metric  graphs  composed  of  n  =  1000  nodes.  For  each  p ,  the  pair 
of  radius  and  energy  providing  minimum  energy  consumption  is 
indicated  by  a  circle  (simulation)  or  a  star  (theory).  We  observe 


c"(r\  _  (9"(r)\2(r)  -  2\'2(r)g  {r))X2(r)  +  g(r)(2\'22 {r)  -  X 2(r)X2(r)) 
{  )  ~  A »(r) 
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Fig.  9.  Second  derivative  of  the  energy  versus  the  transmission  radius  for  sev¬ 
eral  values  of  77. 
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Fig.  10.  Global  energy  consumption  versus  transmission  radius  for  an  RGG; 
theoretical  values  (solid)  and  simulation  (dashed). 


that  the  theoretical  derivations  provide  a  very  good  prediction 
of  the  performance  achieved  by  simulation.  Furthermore,  for 
each  77,  there  is  a  coverage  radius  value  that  minimizes  energy 
consumption. 

The  optimal  values  of  the  coverage  radius  as  a  function  of  77, 
as  predicted  by  our  theoretical  derivations  or  by  simulation,  are 
reported  in  Fig.  11.  From  this  figure,  we  observe  that  there  is 
a  clear  transition  from  the  low  power  attenuation  regime  (i.e., 
77  <  3.5),  where  the  optimal  radius  tends  to  make  the  network 
fully  connected,  as  opposed  to  the  strong  attenuation  situation 
(77  >  4.5),  where  the  network  is  minimally  connected.  This 
means  that  in  the  low  attenuation  case,  minimizing  convergence 
time  is  more  important  than  minimizing  power  consumption. 


i«m 


Fig.  11.  Optimal  transmission  radius  versus  77. 


11  *  \m 


Fig.  12.  Global  energy  consumption  versus  per  node  transmit  power  for  an 
RGG;  theoretical  values  (solid)  and  simulation  (dashed). 


Conversely,  in  the  large  attenuation  case,  it  is  more  important  to 
minimize  power  consumption,  by  limiting  the  number  of  links 
as  much  as  possible,  compatibly  with  the  constraint  of  ensuring 
connectivity. 

Fig.  12  shows  the  average  energy  consumption  versus  the  per 
node  transmit  power  p ,  under  the  same  settings  of  Fig.  10.  The 
circles  (simulation)  and  the  stars  (theory)  represent,  again,  the 
values  of  p  that  minimize  energy  consumption,  for  each  77.  It 
can  verified  the  existence  of  an  optimal  transmit  power  value 
minimizing  the  energy  consumption. 

Finally,  Fig.  13  shows  the  optimal  per  node  power  versus  77. 
For  low  values  of  77,  the  transmit  power  increases  with  77  because 
it  must  cope  with  higher  attenuations  to  guarantee  connectivity. 
Conversely,  for  large  values  of  77,  the  optimal  power  decreases 
with  77  because,  in  such  a  case,  it  is  more  beneficial  to  limit 
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Fig.  13.  Optimal  per  node  transmit  power  versus  77. 


Fig.  14.  Network  algebraic  connectivity  versus  per  node  transmit  power  for 
several  values  of  the  path  loss  coefficient. 


the  number  of  links,  until  the  minimum  power  guaranteeing  the 
connectivity  constraint  is  reached. 

C.  Random  Topologies  versus  Uniformly  Spaced  Grids 

Finally,  it  is  interesting  to  compare  the  energy  consumption 
achievable  with  a  random  geometric  graph  and  a  deterministic 
grid.  In  the  deterministic  case,  we  may  also  distinguish  between 
a  graph  with  nodes  scattered  arbitrarily  in  a  given  area,  but  with 
positions  known,  and  a  regular  uniform  grid,  where  the  points 
are  located  over  a  rectangular  grid. 

1 )  Eigenvalues  of  a  Planar  Uniform  Grid:  We  start  by  de¬ 
riving  the  algebraic  connectivity  of  a  square  grid  whose  n  nodes 
are  uniformly  spaced  within  a  unit  square,  at  a  distance  8  = 
l/m  =  1 1 \fn.  To  avoid  undesired  border  effects,  we  consider 
the  wrapping  of  the  unit  square  in  order  to  form  a  toroidal  sur¬ 
face.  Each  node  is  assumed  to  have  a  link  with  the  neighboring 
nodes  only  if  they  are  at  a  distance  less  than  a  coverage  radius 
r(p),  that  depends  on  the  transmit  power  p,  as  with  geometric 
graphs.  In  Appendix  A,  we  derive  a  closed  form  expression  for 
the  network  degree  and  for  the  algebraic  connectivity  [see,  e.g., 
(45)— (48)]. 

A  numerical  check  of  our  derivations  is  reported  in  Fig.  14, 
where  we  show  the  numerical  value  of  A2  (L) ,  obtained  through 
the  eigendecomposition  of  L,  and  the  value  given  in  (48),  for  dif¬ 
ferent  values  of  77.  From  Fig.  14,  we  can  see  a  perfect  agreement 
between  our  closed  form  expression  and  the  numerical  results. 
Notice,  in  particular,  the  sharp  transition  behavior:  the  eigen¬ 
values  tend  to  n  or  0  depending  on  whether  the  transmit  power 
is  above  or  below  a  threshold;  further,  the  threshold  increases 
with  r\. 

In  Fig.  15,  we  compare  the  algebraic  connectivity  of  the  rect¬ 
angular  grid,  given  by  (48),  with  the  theoretical  value  obtained 
for  the  random  geometric  graph,  given  in  (30),  assuming  the 
same  node  density,  over  the  same  toroidal  surface.  As  shown  in 
previous  works,  see,  e.g.,  [4],  RGGs  tend  to  behave  asymptoti¬ 
cally  as  a  regular  graph.  The  result  shown  in  Fig.  15  is  a  further 
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Fig.  15.  Algebraic  connectivity  of  random  geometric  graph  and  rectangular 
grid  versus  transmission  radius. 


confirmation  of  this  property.  In  the  same  figure,  we  also  re¬ 
port  simulation  results  for  the  algebraic  connectivity,  obtained 
by  averaging  over  100  independent  realizations  of  a  RGG.  We 
can  check,  once  again,  the  good  agreement  between  theory  and 
simulation. 

Finally,  in  Fig.  16  we  compare  the  energy  consumption  ob¬ 
tained  assuming  full  a  priori  knowledge  of  the  nodes’  locations 
or  no  knowledge  at  all.  In  the  first  case,  the  topology  and  the 
power  allocation  over  each  link  are  optimized  according  to  the 
method  illustrated  in  Section  IV.  The  optimal  values,  for  each 
77,  are  indicated  by  colored  dots.  In  the  second  case,  we  report 
(solid  line)  the  energy  consumption  versus  the  average  power, 
assuming  that  all  nodes  transmit  with  the  same  power  (since  they 


35 


n  =  Jtfl 


do  not  have  any  prior  information  about  the  other  nodes  posi¬ 
tions).  To  make  a  fair  comparison  between  the  two  cases,  we 
consider  in  both  cases  the  sum  of  the  powers  associated  with 
each  active  link.  This  means  that,  in  the  second  case,  the  energy 
consumption  is  measured,  for  every  common  transmit  power  p , 

pY^th(p) 

as  £  (p)  =  i=^\  x — ,  where  m  (p)  denotes  the  number  of  neigh¬ 

bors  of  node  i.  The  results  shown  in  Fig.  16  have  been  obtained 
by  averaging  the  results  obtained  in  the  two  settings  over  the 
same  set  of  100  independent  random  deployments  of  n  =  30 
nodes,  uniformly  distributed  over  the  unit  square. 

From  Fig.  16  we  can  see  that,  as  expected,  the  method  as¬ 
suming  perfect  a  priori  knowledge  of  the  node  positions  (in  the 
figure  denoted  as  £(p*))  is  indeed  able  to  achieve  better  per¬ 
formance  than  the  other  method  (£(p)),  as  it  is  able  to  take  ad¬ 
vantage  of  the  a  priori  knowledge.  The  optimal  average  power 
is  the  result  of  a  tradeoff  between  transmit  power  and  conver¬ 
gence  time,  and  always  depends  on  the  path  loss  exponent. 

VI.  Conclusion 

In  this  paper  we  have  addressed  the  problem  of  finding  the 
network  topology  that  minimizes  the  energy  consumption  nec¬ 
essary  to  achieve  consensus  in  a  WSN.  Assuming  a  simple  flat 
fading  propagation  model,  we  have  studied  two  main  network 
models:  1)  arbitrary  networks,  where  the  nodes  are  arbitrarily  lo¬ 
cated,  but  with  known  positions,  and  2)  a  random  model,  where 
the  positions  are  not  known  a  priori,  but  are  modeled  as  random 
variables.  In  the  first  case,  we  have  shown  how  to  optimize  the 
network  topology  and  the  power  allocated  to  each  active  link 
in  order  to  minimize  the  total  energy  necessary  to  achieve  con¬ 
sensus  within  a  prescribed  accuracy.  Topology  optimization  is, 
in  general,  a  combinatorial  problem  and  hence  computationally 
demanding.  To  simplify  the  solution,  we  introduced  a  relaxation 
step  that  enabled  us  to  reformulate  the  energy  minimization 


problem  as  a  convex-concave  fractional  program.  This  alterna¬ 
tive  formulation  can  be  cast  as  an  equivalent  parametric  convex 
problem,  which  enables  efficient  solutions.  The  link  weights  re¬ 
sulting  from  the  solution  of  the  convex  problem  must  then  be 
thresholded  to  find  out  the  network  topology.  We  have  shown 
through  numerical  results  that  the  thresholding  operation  may 
yield  a  considerable  reduction  of  the  number  of  active  links,  yet 
with  very  limited  effect  on  performance. 

In  our  deterministic  setting,  we  assumed  point-to-point  links 
and  we  optimized  the  power  over  each  link.  In  practice,  a  con¬ 
sensus  algorithm  running  over  a  wireless  network  could  ben¬ 
efit  from  the  broadcast  channel.  Deciding  between  broadcast  or 
one-to-one  links  entails  a  proper  choice  of  the  medium  access 
strategy,  to  establish  in  which  time  slot  each  node  has  to  listen 
to  which  broadcaster. 

Conversely,  in  the  random  network  case,  lacking  any  infor¬ 
mation  about  the  internode  distances,  we  assumed  a  broadcast 
communication  strategy.  In  such  a  case,  the  network  topology 
is  modeled  as  a  random  geometric  graph.  We  have  derived 
closed  form  expressions,  albeit  valid  only  asymptotically,  for 
the  algebraic  connectivity,  assuming  a  common  transmit  power. 
Then,  building  on  these  expressions,  we  have  shown  that  the 
energy  consumption  is  a  convex  function  of  the  coverage  radius. 
We  have  also  shown  that  a  random  geometric  graph  performs, 
asymptotically,  as  a  regular  graph  built  over  a  rectangular  grid. 
This  confirms  previous  results,  although  now  in  the  context  of 
energy  minimization  over  consensus  networks. 

Finally,  we  have  compared  the  performance  achievable  with 
arbitrary  and  random  graphs.  Clearly,  the  knowledge  of  the 
node  locations  allows  better  power  allocation,  that  translates 
into  lower  energy  consumption  to  achieve  consensus.  However, 
in  practice  there  is  a  price  associated  with  the  knowledge  of 
node  location.  This  knowledge  requires  the  acquisition  of 
the  node  positions  first  and  then  a  centralized  optimization. 
Conversely,  the  random  approach  can  be  followed  also  in  a 
decentralized  fashion,  with  only  minimal  information  about 
some  global  parameters  like  number  of  nodes  and  area  covered 
by  the  network,  and  it  does  not  need  any  extra  hardware  or 
computation  to  acquire  the  nodes’  locations. 

In  this  paper,  we  have  assumed  a  simple  flat  fading  channel 
model,  whose  effect  is  only  to  introduce  attenuation  and  su¬ 
perimpose  noise.  The  simple  model  captures  the  essence  of  the 
problem  and  keeps  the  overall  problem  complexity  under  con¬ 
trol.  However,  looking  at  potential  applications,  it  would  be 
interesting  to  generalize  the  approach  to  the  case  where  the 
channel  model  is  more  complicated.  Furthermore,  the  network 
topology  has  been  assumed  to  be  static.  However,  a  dynamic 
topology,  possibly  adapted  to  the  consensus  state,  might  pro¬ 
vide  better  performance. 


Appendix  A 

In  this  section,  we  derive  an  analytical  expression  for  the 
eigenvalues  of  the  adjacency  matrix  A  of  a  rectangular  grid  over 
a  unit  torus,  as  a  function  of  the  coverage  radius  r(p).  We  as¬ 
sume,  for  simplicity  that  the  number  of  nodes  n  is  a  square 
number,  i.e.,  n  =  n\,  with  n\  integer.  We  number  the  rows  of 
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the  grid,  1  to  ni,  going  from  bottom  to  top.  The  n  x  n  matrix  A 
can  be  expressed  in  a  cyclic  block  form,  as 


^4-0 

Ai 

a2  . 

•  •  Ani  _2 

Ani  —  1 

Ani  —  i 

Ao 

Ai  . 

•  •  —3 

An  i_2 

A  = 

^-ni-2 

A-n1  —  1 

Ao  . 

•  •  4Lni  _4 

Aii  —3 

- 

a2 

A3  . 

•  •  4l_ni  _  1 

^4-0 

where  each  block  Ai,  for  i  =  0, ,  * . ,  n\  —  1,  is  an  n\  x  n i 
circulant  matrix.  In  particular,  the  block  Ao  is  the  adjacency 
matrix  between  the  nodes  lying  on  the  same  row  of  the  grid, 
starting  from  the  leftmost  node  and  proceeding  along  the  row. 
Ai  represents  the  adjacency  matrix  between  the  nodes  of  row  j 
and  the  nodes  lying  on  the  (i  +  j)th  row,  for  i  =  1, . . . ,  ni  —  j, 
and  on  the  {i  +  j  —  ni)th  row  for  i  =  n\  —  j  +  1, , . . ,  n\  —  1. 
The  toroidal  structure,  used  to  avoid  border  effects,  reflects  into 
the  cyclic  structure  of  A,  so  that  we  have  Anx~i  =  Ai  if 

i  =  (m  —  l)/2  +  1, . . . ,  ni  —  1  and  ni  odd 
i  =  ni/2  +  1,  < . , ,  m  —  1  and  even. 

We  also  note  that  A  is  a  block  circulant  matrix  with  circulant 
blocks  (BCCB).  Since  each  block  Ai  for  i  =  0, . . . ,  n\  —  1  is 
anni  x  n\ -dimensional  circulant  matrix,  the  eigenvalues  of  Ai 
admit  the  following  decomposition 

Ai  =  FniAiFni 

with 

1  ' 

u ni~1 


(j 

where  u  =  exp  (—  jlit/rii)]  A i  is  a  diagonal  matrix  whose  en¬ 
tries  are  the  eigenvalues  of  Ai,  i.e.,  the  discrete  Fourier  trans¬ 
form  (DFT)  of  the  first  row  of  Ai.  Hence,  the  matrix  A  admits 
the  following  diagonalization: 


F  — 

1  ni  — 


V^T 


l  l 

1  (J 


1  a/11-1  uoni~2 


for  /,  m  =  0, . . . ,  ni  —  1  and  the  coefficient  (p)  is  £;th  entry 
of  the  first  row  of  Ai  for  fc,  i  =  0, . . . ,  n\  —  1.  Equivalently,  (37) 
can  be  seen  as  the  two-dimensional  DFT  of  a  matrix  A  whose 
rows  A(z, :)  are  given  by  A(i, :)  =  A^(l, :)  for  i  =  0, . . . ,  n\  — 
1,  applying  the  1-DFT  to  the  rows  and  then  to  the  columns  of  A. 

In  general,  for  a  regular  grid  where  each  node  has  degree  d(p) , 
the  algebraic  connectivity  can  be  expressed  as 

A2  {L{p))  =  d(r(p )) 

n±  —  1  ni  —  1 

-  E  E  a£\p)exp(-j2Trk/n1)  (38) 

i=0  k= 0 

where  the  coefficients  a^\p)  can  be  expressed  as  follows.  Let 
u(x)  denote  the  step  function,  i.e.,  u{x)  =  1  for  x  >  0  and  zero 
otherwise.  For  m  odd,  we  have 

aQ°^(p)  =  0  and 

ao  \p)  =  u(r(p)  -  i6)  i  =  1, . . . ,  («1  -  l)/2,  (39) 

(p)  =  u(r(p)  —  kS)  and 

ak\p)  =  u(r(P )  _  \/*2  +  k?6)  i,k  =  1, . . . ,  («i  —  l)/2 

(40) 

ak\p)  =  aS-k(P)  k=(m-  l)/2  +  1, •  •  • , «1  -  1  and 

*  =  o,...,(m-i)/2.  (41) 

Finally  Ai  =  Ani-i  for  i  =  (n i  —  l)/2  +  1, . . . ,  ni  —  1. 
Similarly,  in  the  case  of  ni  even,  we  have 

40)(p)  =  0  and 

ao\p)  =  u(r(p )  ~iS)  i  =  !)•••) ™i/2>  (42) 

a  I'0-1  (p)  =  u{r{p)  —  kS)  and 

ak\p)  =  u(r(p)  ~  "\A2  +  fc2  8)  i,  k  =  1, , . .  ,ni/2,  (43) 

4°  (P)  =  an\-k(P )  k  =  ™i/2  +  1,  •  •  • ,  ni  -  1  and 

*  =  0,...,m/2  (44) 

with  Aj  =  A^-;  if  *  =  rti/2  +  1, . . . ,  -  1. 

Finally,  we  can  derive  analytical  expressions  of  the  network 
degree  and  of  the  algebraic  connectivity  for  m  odd  as 


A  =  (Fni  <8> Fni)HA(Fni  ®Fni) 
where  (8)  denotes  the  Kronecker  product  and 

ni  —  1 

A=  ^  fc 

fc=0 

with  fini  =  diag[l,  a;, . . . ,  a;ni-1].  The  entries  A/?m(A)  of  the 
matrix  A  can  be  also  expressed,  after  a  few  algebraic  manipu¬ 
lations,  as 


n\  —  1  n  i  —  1 

Vm(A)  =  E  E  4°(p) 

2=0  fc=0 

x  exp  (—  j2'Kklfn-])  exp  (—j27rmi/ni)  (37) 
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and  for  n,\  even 
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x  (1  -  cos(27rfc/rai))  +  4  53  ani/2(p) 

2=1 

+  2ai01)/2(p)  +  2o£/22)(p) 
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2  =  1  fc=l 

For  example,  in  the  case  where  r(p)  <  \[2b,  i.e.,  assuming  a 
regular  grid  of  degree  d(p)  =  4,  we  obtain 


Ani  _i  =  Ai  =  I  and  Ai  =  O  Vz  =  2, . . .  ,ni  —  2  (49) 

while  Ao  is  a  circulant  matrix  whose  first  row  is  the  ni -di¬ 
mensional  vector  [0100 ...  001].  Then,  the  eigenvalues 
are  given  by  A/(A0)  =  2 cos(27r//ni)  and  applying  the 

formula  in  (37)  the  eigenvalues  of  A  are  A/?m(A)  = 

2  cos(27r//ni)  +  2  cos(27rm/ni)  for /,  m  =  0, . . . ,  ni  —  1. 
The  second  largest  eigenvalue  of  A  is  obtained  for  l  =  1  and 
m  =  0  and  is  given  by  Ai?q(A)  =  2cos(27r/ni)  +  2.  Thus 
the  eigenvalues  of  the  Laplacian  matrix  can  be  expressed  as 
A i(L)  =  4  —  A; (A)  and  the  algebraic  connectivity  of  the 
network  is  given  by  A2  (L)  =  2  —  2  cos(27r/ni). 
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